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There is considerable evidence for some form of charge 
ordering on the hole-doped stripes in the cuprates, mainly as- 
sociated with the low-temperature tetragonal phase, but with 
some evidence for either charge density waves or a flux phase, 
which is a form of dynamic charge-density wave. These three 
states form a pseudospin triplet, demonstrating a close con- 
nection with the E®e dynamic Jahn- Teller effect, suggesting 
that the cuprates constitute a form of Berryonic matter. This 
in turn suggests a new model for the dynamic Jahn- Teller ef- 
fect as a form of flux phase. A simple model of the Cu-O bond 
stretching phonons allows an estimate of electron-phonon cou- 
pling for these modes, explaining why the half breathing mode 
softens so much more than the full oxygen breathing mode. 
The anomalous properties of 2 ~ provide a coupling (corre- 
lated hopping) which acts to stabilize density wave phases. 



I. INTRODUCTION 

In a (molecular) Jahn- Teller (JT) effect, an electronic 
degeneracy in a symmetric molecule is lifted by distorting 
the molecule and lowering the free energy by an amount 
Ejt- A dynamic JT (dJT) effect arises when there is 
a degeneracy both of electronic levels and of molecu- 
lar vibrations (strictly speaking, when two independent 
distortions have the same JT energy lowering). Quan- 
tum tunneling between the two distortions restores the 
(time-averaged) symmetry of the molecule, but leads to 
a net orbital angular momentum which is half-integrally 
quantizedElu. In principle—an analog of the dJT effect 
should exist in dense solidsma, but this feature of ground 
state nuclear motion has not been adequately described. 
Recently, a number of proposals have been made for solid 
state systems with finite orbital angular momentum in 
their grouttdstates, including 'Berryonic roatter'El and the 
flux phaseQ u or orbital antiferromagnetliJiii Here the 
connection between these concepts is worked out, a for- 
malism is presented for describing the dJT effect in terms 
of phonon-assisted flux phases. 

This result is particularly timely, since dJT effects have 
been proposed to play a role in a number of materialS|-af 
current interest ..including buckyballstL3i3, cupratesEij E.3, 
and manganitesEEI. Remarkably, many of these materials 
display stripe phases and/or high temperature supercon- 
ductivity, and it has bee^-proposed that either strong 
electron-phonon couplingcHtJ or flux phases may be re- 
sponsible for these features. Moreover, in the cuprates, 



there is recent evidence for . bot h r-the flux phase@ and 
for strong phonon anomaliest§c2rc3. Hence, the present 
formalism is applied to a simple model of the in-plane 
phonons in the cuprates. 

The two JT modes plus the dJT state can be combined 
into a 'pseudospin' tripletc3 of site charge-density wave 
(CDW), bond CDW, and flux phase (dynamical CDW). 
Restricted to a single four-copper plaquette 'molecule', 
they reduce to an E® (61 + 62) JT problem. The site and 
bond CDW's can be related to important phonon modes 
in the cuprates, namely the oxygen breathing mode and a 
shear mode which couples to the tilting mode associated 
with the low-temperature tetragonal (LTT) phase. Since 
there is considerable interest in the flux phase, some care 
is taken to analyze it in terms of the associated struc- 
tural distortions. In particular, the anomalous n phase 
is just the Berry phase of the dJT effect. In the molecu- 
lar problem, the dJT effect arises at a (weak) degeneracy 
point; remarkably, at this point the dynamics is (classi- 
cally) chaotic. 

Our model calculations find a crossover from an anti- 
ferromagnetic insulator at half filling to a doped para- 
magnetic phase stabilized by density wave order. (The 
detailed doping dependence is probably complicated by 
the presence of stripes, as discussed in the Appendix.) 
The density wave is predominantly stabilized by electron- 
electron interaction, but with a significant structural 
component. The system is close to but probably not at 
the dJT degeneracy point, with the shear and flux phases 
lying lowest in energy. 

This paper is organized as follows. Section II surveys 
the various phonons which experimentally are found to 
couple to the holes, and shows that they fall into three 
classes, including site and bond CDW's and the flux 
phase. It is explained how the bond CDW couples to 
LTT and dimpling distortions. (In an Appendix addi- 
tional evidence is presented that the coupling is partic- 
ularly strong in the hole-doped regions, i.e., on charged 
stripes.) Section III shows that these classes have a well- 
defined group (pseudospin) structure (Section IIIA), and 
are closely associated with the dJT effect of a square 
molecule (IIIB). Section IIIC generalizes the pseudospin 
formalism to a lattice. In Section IV these results are ap- 
plied to the cuprates, using both a three band (IVA) and 
a one band (IVB) model to estimate the electron-phonon 
coupling of a large number of planar bond stretching 
modes, including some outside of the above three classes. 
The reason why the half breathing mode is so strongly 
renormalized is explained. A simple model of phononic 
flux phases is presented which can provide a basis for 
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a general theory of band dJT effects. In Section V 
the self-consistent gap equations are solved for combined 
electron-electron and electron-phonon coupling, and it is 
found that the anomalous behavior of the 2 ~ molecule 
leads to a particular form of correlated hopping which 
helps stabilize paramagnetic CDW phases. Some discus- 
sion of desiderata for a more complete model are pre- 
sented in Section VI, and conclusions are presented in 
Section VII. 



II. CLASSIFICATION OF PHONON MODES 

There is considerable evidence for strong electron- 
phonon coupling effects in the cuprates. While the ef- 
fects are weaker than in the nickelates and manganites, 
we postulate that it is this relative weakness which al- 
lows superconductivity to compete against charge order- 
ing. The difficulty in analyzing the coupling lies in the 
fact that it is spread over so many different modes, and 
the modes vary from cuprate to cuprate. Here we show 
that the strongly coupling modes can be classified into 
three groups: site and bond CDW's and flux phases. 
Additional evidence specifically tying the phonons to the 
charged stripes is presented in the Appendix. 

The three categories of phonons can be described in dif- 
ferent ways. Perhaps the simplest is as site CDW (lead- 
ing to pileup of charge on alternate Cu's) vs bond CDW's 
(pilcup on oxygens) vs dynamic interchange of the two 
(the flux or orbital antiferromagnetic - OAF - phase). 
Alternatively, since the modes all couple strongly to a 
Van Hove singularity (VHS), the first two may be clas- 
sified by how they split the VHS peak: the site CDW's 
leave the X and Y ((tt, 0) and (0, tt)) point VHS's equiv- 
alent, but split each (as in the one-dimensional Peierls 
transition); the bond CDW's make the X and Y points 
inequivalent. 

A. Breathing Modes 

Strong coupling to oxygen breathing modes is found 
in a number of cuprates. WeberEZI originally suggested 
that there should be large coupling of the electrons to 
the oxygen breathing mode at (jr, tt). However, this 
mode is suppressed near half filling by strong coupling 
effects, since the CDW requires a charge imbalance on 
different copper atoms. The coupling is enhanced by 
proximity to the VHS, and hence-muldJ.ead to a CDW 
instability near optimal dopingoEZrcj'Ej. Experimen- 
tally, breathing modes appear more strongly coupled near 
(tt, 0) than at (n,n), both-in La2_ :r Sr. r Cu04 (jLSCO) 
and YBa 2 Cu 3 07-5 (YBCO)EI. McQueeney, et alE9 have 
found that the phonon spectra show evidence for a dy- 
namic period doubling in the T — (tt, 0) direction, which 
they associate with stripes. Since this period is indepen- 
dent of dopingcil, it is unlikely to be related to stripe 



periodicity, but could be associated with charge ordering 
along the stripes. A possible explanation for the strong 
coupling of this half breathing mode is given in Section 
IVA. Very similar anomalies are found in the nickelates, 
with phonon softening and splitting over a large part of 
the Brilloiiin zone, including evidence for local modes tied 
to stripesEi In LSCO, stripe correlations have a sfmng 
effect on these oxygen-related phonons near 70 meVcdS. 

In this same energy range (50-70meV_below the Fermi 
level), photoemission measurementsoEj have found a 
break in the dispersion of the cuprates, present through- 
out A;-space, although having much more dramatic effects 
near (tt, 0). While there-are alternative interpretations^, 
it has been proposedlla that this is the conventional 
renormaliaation associated with strong electron-phonon 
couplingo to the breathing mode phonons near this en- 
ergy. Such renormalization can .produce the extended 
saddle points observed near (tt, 0)eB. 

B. Relation of Shear Mode to LTT, Dimpling, etc. 

Evidence for proximity to a shear mode instability has 
been noted in LSCO near the 1 /8th anomalyy-in the form 
of softening of the C\\ — C\2 elastic modulus". However, 
it is remarkable that, whereas theoretically there should 
be a strong electron-phonon coupling with bond stretch- 
ing modes, the observed soft modes are bond bending 
modes, with different modes in each family of cuprates: 
the low-temperature orthorhombic (LTO) and LTT tilts 
in LSCO, a dimpling mode in YBCO, and a rotation 
of the four planar oxygens about the central Cu in the 
electron-doped cupratesEZl. Here, we suggest that there 
is a connection between these modes and the shear insta- 
bility. 

We first consider LSCO. It is a common finding that 
the local structure in the cuprates differs from the av- 
erage structural order determined by neutron or x-ray 
diffractionta, and it has been suggested that the observed 
LTO phase in LSCO is really dynamically disordered, 
with the local order being closer to LTTEJl. [There have 
been numerous hints of local LTT orderEJCJ, but the 
problem is complicated by the simultaneous occurence of 
stripe order. The issue of the relation between stripe or- 
der and the LTT distortion is addressed in the Appendix. 
Here, we merely note that the relevant local distortion is 
likely to be of LTT form, and explore the connection of 
this distortion to the shear mode. 

By increasing the Cu-0 bond length along the tilted 
bonds, both the shear and the LTT distortions make the 
two oxygens in a unit cell crystallographically inequiva- 
lent, thereby coupling to the VHS. Such, VHS splitting 
has recently received renewed interestc3'c3. We have sug- 
gested that the LTT tilt is a secondary modification oi-a 
shear mode instability due to lattice mismatch effectsEJ: 
the underlying instability is to a shear (making inequiv- 
alent bondlcngths), but the Cu02 planes are under com- 
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pression, so the LTT tilt arises as a compromise, making 
the Cu-0 bondlengths inequivalent, without increasing 
the net Cu-Cu bondlength. Qualitatively, the LTT tilt 
pattern arises because the CuC>6 octahedra tilt as rigid 
unitsell: since the octahedra are corner sharing, alternate 
octahedra must tilt in opposite directions, leading to a 
distortion at wave vector Q — (tt,tt). (The detailed in- 
terpretation involves partial hole localization on oxygens, 
and the large size difference between 0~ and 2 ~£3) 

Similar effects arise in other cuprates. Thus, high 
energy phonons behave similarly in LSCO and YBCO, 
but at low frequencies there are very different struc- 
tural distortion. In YBCO there is little evidence for 
an LTT phase (although there was a report of a local 
tilt- mode instability^). More likely, the role of the LTT 
distortion is taken by the dimpling mode. While this 
mode involves predominantly copper and oxygen mo- 
tion out of the Cu02 planes, it also couples to in-plane 
Cu-0 displacements. In particular, it couples to the 
orthorhombicityEJ, and hence to the same shear mode 
involved in the LTT phase (distortion Bi in Fig. [2] be- 
low). Hence theilimpling also splits the VHS's, but in a 
different patternEj. Experimentally, the distortion com- 
petes with superconductivity, and a large decrease of the 
orthorhombicity is found both at and below TJE3, while 
a discontinuous change in the dimpling is found on pass- 
ing through optimal dopingpl. Chain ordering also plays 
a role, since it splits the degeneracy of the VHS's along 
and across the chain direction. Photocmission seems to 
find that the VHS along X (Y) is below (above) the 
Fermi levels, while neutron scatteringtil finds that the 
stripes align along the chain direction, reminiscent of 
LTT pinning. In a related 123 compound, (Lai-aCa^)- 
(Bai.75_a;Lao.25+a:)Cu30y, the La 3+ at the Ba cite dis- 
rupts the chains, leading to an average tetragonal struc- 
ture; in this case a plane buckling is observed, which is 
largest at the doping of optimal superconducting T c xB 

A similar effect is found in Bi2Sr2CaCu20s+,5 (Bi2212). 
The 'ghost' Fermi surfaces in Bi2212 are clear evidence 
for a lattice period doubling, and a connection with the 
antiferromagnetic order has been suggested^. However, 
since they are present at room temperature and at op- 
timal doping - far from the region where strong mag- 
netic correlations would be expected - they may instead 
be related to a structural transition - the known or- 
thorhombic orderingEll. This lattice anomaly also couples 
strongly to the VHS. Most suggestively, it has recently 
been found that the intensity of these ghost features peaks 
near optimal doping^. This is close to the prediction 
of the VHS stripe model, in which the structural order 
should be maximal in the slightly overdoped regime, just 
when magnetic correlations disappear (termination of the 
stripe phase). 



C. The Flux Phase 

The flux phase has generated intense interest because 
of its highly unusual physical properties. The experi- 
mental evidence for this phase, which remains contro- 
versial, includes (i) the fit of the dispersion of the Mott. 
insulator Sr 2 Cu0 2 Cl 2 (SCOC)El to the flux phase formtl 
(however, alternative interpretations include a Mott anti- 
ferromagnet, while jaorc recent experiments find a more 
complex dispersionEj), £ii) the d- wave-like nature of the 
normal-state pseudogapEH, which has been interpreted in 
terms of the conical point (node) in the flux phase dis- 
persion at (7r/2, 7r/2), and (hi) the recent observation of 
the induced magnetic fields associated with the pxbital 
currentalj. (However, more recent measurementsEj find 
that the magnetic moment lies in-plane, and not along 
the c-axis as originally expected, and the experiment also 
has alternate interpretations^!) 

D. Competition and Cooperation of the CDW's 

In Section IIIB, we will show that the molecular equiv- 
alent of the flux phase, the dJT state, arises from com- 
petition between the two other modes. In the cuprates, 
there is also evidence for a competition between the two 
types of CDW, which we briefly review here. In LSCO, 
both the breathing and the LTT modes show anomalous 
softening with doping. An interesting situation arises in 
YBCO. Pintschovius and Reichardt (p. 349 of Ref. |(]) 
note that hole doping may renormalize the half-breathing 
mode phonon frequency from 74meV_to ~37meV, due to 
a special form of Van Hove nestingpU. The mechanism 
proposed involves VHS induced coupling of two phonon 
modes remarkably similar to the flux phase effects we 
are describing. The dimpling of the Cu02 planes causes 
the VHS's to bifurcate (e.g., the X VHS splits along the 
r — X — r line). In turn, nesting between the bifurcated 
peaks leads to breathing mode softening along T — X 
and can lead to a period doubling instability. The ex- 
act nesting wave vector is controlled by the degree of 
VHS bifurcation, which in turn depends on the amount 
of dimpling. Thus, two phonon modes (here the dim- 
pling and half-breathing modes) are strongly coupled by 
the VHS. A similar competition between breathing and 
octahedpal tilt modes near a superconducting instability 
is foundEj in Bai-^K^BiOs; in this case the tilt modes 
are near the R point of the Brillouin zone, coupling to 
the three-dimensional VHS's. 

Other routes to mode coupling are alsojiossible. For 
example, when the breathing mode softensEj'EjO, it can 
couple to a lower frequency phonon, either an LTT distor- 
tion in LSCO or a dimpling mode in YBCO, and drive 
that second mode soft. This strong mixing of the two 
modes can create a flux phase. 



3 



III. ELECTRON-PHONON COUPLING: FROM 
SQUARE MOLECULES TO CUPRATES 



Brillouin zone. This .same problem arises in SO(5)E3, and 
the Henley operatora was introduced to deal with it. 



In this section, we show that the three phonon modes 
form a symmetry group, and map the interaction to a 
molecular JT model. We demonstrate that the same de- 
generacy arises in both the molecular and lattice prob- 
lems. In Section IV we apply these results to the 
cuprates, showing how the phonons in LSCO fit into the 
scheme, and calculating their coupling to electrons. 



(a) (b) (c) 




FIG. 1. Pseudospin triplet: Ocdw (a), r (b), Ooaf (c). 



A. The Pseudospin Group 

We begin with a purely theoretical description, based 
on the gpoup theory of instability, or spectrum generating 
algebrasE3, specialized to the VHS instability grouptilLj, 
of SO (8) symmetry. Within this group there is a triplet 
of CDW operators, which couple strongly to phonons. 
These operators are spanned by an SU(2) subalgebra 
which we call 'pseudospin'. 

To specify the order parameters in terms of bilinear 
electronic operators, it is necessary to define these op- 
erators on a larger unit cell, containing a plaquette of 
four copper atoms. It is convenient to introduce creation 
operators cJ ± a (r) = {i>^ i !T (r)±i>^ IJ {r))/2. The pseu- 
dospin can be written in terms of the Fourier transforms 
of these states. In the basis 

{B a (k)} = {a +>a (k),a^(k)}, (1) 

the explicit form of these matrices is 

*=\ f *, ( 2 ) 

Ooaf = ^ T y , (3) 

Ocdw = ? z , (4) 

where 

*.-(; o). < 6 > 



The electronic arrangements corresponding to these 
operators are illustrated in Fig. [l]: Oc dw is a site CDW, 
with excess charge on the odd sites 1 and 3, and a deficit 
on the even sites 2 and 4; r is a bond centered CDW, with 
excess charge on the 1-2 and 3-4 links, and a deficit on 
the 1-4 and 2-3 bonds; while Ooaf is an orbital antifer- 
romagnet, with an orbital current flowing clockwise, from 
1 ^ 4 ^ 3 ^ 2. In a three-band model, Ocdw would 
represent a Cu-based CDW, and r an O-based CDW. 
Note that the CDW mode is equivalent tOpthe charge 
ordering found in La2™Sr x Ni04 for x ~ 0.5E2. The or- 
bital current operatoiO Ooaf is closely related to the 
flux phaseH. In the flux phase, hopping is accompanied 
by an accumulation of a phase change by n on circulating 
a plaquette. For the plaquette of Fig. Qc, this leads to a 
contribution to the Hamiltonian proportional to 

H' = e"' '^{a\ai + a\a,2 + a\a% + a\a,4) + h.c. 

= —j={a\ai + a\d2 + a\a^ + a\a,4 + h.c.) — 2\f2T y . (8) 
V2 

The first term in Eq. || corresponds to an (uninteresting) 
uniform hopping on the plaquette. Thus, up to this term, 
the flux phase is equivalent to the orbital antiferromagnet 
Ooaf (T y ). 




with Ci = cos kia and 7 = c x — c y . The operators Ocdw, 
f , and Ooaf transform exactly as a pseudospin near the 
VHS's (c x — c y = ±2), but not at a general point in the 



FIG. 2. Bi (a) and B2 (b) distortions of CU4O4. Open 
(filled) circles represent the Cu's (O's). 

While the instabilities are predominantly electronic, 
they are accompanied by lattice instabilities, as in the 
Peierls transition. In particular, the CDW couples to 
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the oxygen breathing mode and the r-mode couples nat- 
urally to a shear distortion; these are referenced to a 
single molecule in Fig. ||. The relation of this shear to 
the LTT tilting mode was clarified in Section II. 

This pseudospin triplet is an exact condensed-matter 
analog to the molecular JT modes responsible for the dJT 
effects found in molecules. The enlarged unit cell asso- 
ciated with Eq. [I] provides a natural connection. Thus, 
we regard the CuC>2 planes as composed of a square ar- 
ray of CU4O4 'molecules' bonded to each neighbor by 
a pair of oxygens. These molecules can be embedded 
into the lattice in two equivalent ways, as illustrated in 
Fig. [| While it is simpler to analyze a square lattice of 
square molecules, as in Fig. ||a, a more accurate repre- 
sentation is in terms of a 'centered square' lattice as in 
Fig. |^b. [Crystallographically, the centered square is not 
a primitive unit cell, but it is convenient for the present 
purposes.] The latter case has twice as many molecules, 
but they are corner sharing. The atoms (coppers and 
oxygens) in each molecule deform in the same way. For 
instance, Fig. ||c illustrates the flux phase, labelling each 
square plaquette which contains a +7r flux with a plus 
sign. When interpreted as in Fig. ||b, the ground states 
have a two-fold degeneracy, which is also found in the 
flux and CDW phases. 



(a) 

nn 
nn 



(b) 



(c) 



F 




R 




+ 











FIG. 3. Berryonic matter: two views of the lattice of square 
molecules, either as simple square (a) or as centered square 
with corner sharing molecules (b) compared with the flux 
phase (c), where the +'s indicate plaquettes with a +7r flux. 

The above analysis resolves an old conundrum. Earlier 
calculations interpreted the LTT phasedn terms of a Van 
Hove - Jahn-Teller (VHJT) effect£§E3, where the elec- 
tronic degeneracy is associated with the presence of two 
VHS's. This is puzzling, since JT effects usually have a 
simple molecular basis, but such a basis is not obvious in 
the present case. We now see that the VHJT effect on 
the conventional lattice is equivalent to a conventional 
JT effect on a supercell lattice. The resulting lattice con- 
stitutes a nearly exact realization of a Berryonic crystal 
- replacing the originally proposed!] JT active triatomic 
molecules by square molecules. 

The close connection between the JT and Van Hove 
viewpoints can be clarified by representing the CuOa 
plane electronic states in a supercell representatio n! 1 ^P 6 ! 
replacing the copper (or Zhang-Rice singlet) wavefunc- 
tions by electronic states symmetrized on a single plaque- 
tte. These states ip+ijk — (ipi + ^2 +iV'3 + kip^/2 have 
Mg (i>++++), B 2g 00+-+-), and E u — , tp+ — +) 
symmetry. Including nearest (t) and next-nearest neigh- 



bor (t') hopping, the kinetic energies are 

E{A lg ) = -2t-t', 
E{E U ) = t', 
E(B 2g ) =2t- t', 



(9) 



so at half filling, the Ai g level is filled and there are two 
electrons in the E u levels. Optimal doping corresponds 
to one extra hole per plaquette, leaving one electron in 
the E u levels, and a JT effect. The full band disper- 
sion of the square lattice can be written as a superposi- 
tion of these four states (Appendix I of Ref. |l5|). It is 
found that all states near the VHS's are built up exclu- 
sively of E u states. On the other hand, states near the 
nodal points have A\ g character. Note that this explic- 
itly demonstrates that VHS instabilities persist down to 
'lattices' as small as a single 2x2 plaquette. 

The above calculations are for the weak coupling limit, 
and are consistent with the charge ordered states found 
in Section V. However, the analysis can be repeated in 
the strong coupling limit. Thus, for an antiferromagnetic 
arrangement, the (mean field) energies are E{A\ g ) 

E(B 2g ) 



f - VW r T4£{ T -E{E u 



t' - U, with 



U ~ U/2. With parameters^ t ~ 325meV, U/t ~ 6, 
t'/t ~ -0.276, and J = 4t 2 /U, E(E U ) - E(A lg ) ~ 
2t' + J ~ t/6 > 0, so the first hole would again come 
from the E u level, although all levels are close in energy. 
(Note that a smaller J would reverse the order of the 
levels but would still lead to an electronic degeneracy.) 
Thus, one can still have a JT instability in the strong 
coupling limit, although its nature can be significantly 
modified. For a somewhat different example, see Ref. 7CL 
The present JT effect is significantly different from the 
conventional JT effect of a CuOg octahedron, which in- 
volves the electronic (pseudo-) degeneracy associated with 
the Cu d z 2 and d x 2_ y z orbitals. Many earlier studies 
(e.g., Refs. [T^, [n]) have proposed that this (pseudo) JT 
effect plays an important role in cuprate physics. The 
present model works even if the d z 2 orbitals are com- 
pletely uncoupled (although a residual coupling could 
enhance the presen t effe cts). The model is closer to the 
large polaron limitOO, but with a particular choice of 
polaronic coupling dictated by the underlying SO(8) sym- 
metry. 



B. The Molecular Jahn-Teller Effect 

As a preliminary to studying the VHJT effect on a 
lattice, we first look at the molecular analog, a CU4O4 
molecule with square planar symmetry, D^^Mhich cor- 
responds to the E <g> (61 + 6 2 ) JT problerxfl'Bfca. This 
molecule can display a dJT effect which bears a strik- 
ing resemblance to the flux phase. Figure || displays 
the relevant B\ and B 2 distortions, assuming predom- 
inantly oxygen vibrations. The B2 distortion is the oxy- 
gen breathing mode, while the Bi has the form of the 
shear wave associated with the LTT instability. 
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The JT energy can be written E% = V?/(2Mwf), 
where Vi is the electron-phonon coupling and M is the 
ionic mass. In terms of Ejt, there are three classes of 
solution: 

(i) E JT ^ E JT-> 



(ii) E% = E$ but ui ^ 



p(2) 



(iii) Ejj, — Ej Z ^ and lj\ = u>2- 
For case (i) the lowest energy state consists of a static 
JT distortion of the mode with larger JT energy only. 
This simple case may apply to the cuprates: it has long 
been a puzzle why the oxygen breathing modes do not 
display the large softening expectedEH near (n, n). On the 
other hand, the LTT mode is quite soft, and nearly un- 
stable. This would follow if E$ > Ef T . The results of 
Section V are consistent with this possibility. However, 
formally the special cases (ii) and (iii) are more interest- 
ing, allowing dynamic solutions which are closely related 
to the flux phase. 

Case (iii) reduces exactly to the,well-known E®e prob- 
lem of the triatomic molecule»EIO, and the electronic 
wave function is double valued: when the molecular mo- 
tion undergoes a 2?t rotation, the electronic wave 
tion picks up an extra factor of it, the Berry phase 
However, in the cuprates the phonon frequencies are very 
different (Section IV), so case (iii) does not apply. Nev- 
ertheless, condition (ii) (weak degeneracy) may hold ap- 
proximately in the cuprates, and it is known that inclu- 
sion of quadratic vibrpaic coupling enhances the range 
of case-(ii) deeeperacyll3. This special case has recently 
been analyzedO: the dJT effect with n Berry phase is 
preserved, although the classical motion may be chaotic. 



,(2 



m. 




FIG. 4. Dynamic JT rotation mixing B± and B2 modes 
of O4. Large open circle represents one hole. Each time step 
(a-f ) represents one quarter of a cycle of molecular rotation, so 
in frame e the molecules have completed one cycle of rotation 
with respect to frame a, while the hole has completed only 
half a rotation. 

The origin of the wave function sign change is ex- 
plained in Fig. ||. As each oxygen atom rotates about its 
undistorted position, the hole (large circle; a second hole 
is in the symmetric position) rotates about a different ori- 
gin, the center of the square. When the individual atoms 



complete one full rotation, the hole has only completed 
half of a rotation, just as in the triangular moleculea. In 
analogy with the triatomic molecule, the electronic wave 
function can be kept single valued by multiplying it by a 
phase factor, e l9 / 2 . 

Thus, the dJT effect with anomalous Berry phase ex- 
ists for a square molecule; while some 'fine tuning' of the 
parameters is required, this may arise from quadratic vi- 
bronic couplingjla. In the following subsection, we will see 
that a similar accidental degeneracy exists for the lattice 
problem, and it is known that rthe-structural transitions 
in LSCO are strongly nonlineartilCa. However, in Section 
VC it will be seen that on a lattice the flux phase can 
be stabilized over a wider parameter range, even without 
nonlinear coupling. 



C. Extension to the Lattice 

Thc_E (g) (61 -I- 62) problem has been extended to the 
latticeEJ but only the nondegenerate case (i) has been 
analyzed in detail. Here, we demonstrate that the de- 
generate case (ii) remains a possibility. On a lattice, the 
JT coupling becomes 



H 



JT 



»=1,2 



H 



JTA 



(10) 



T q m 

(11) 

where the m-sum is over different phonon branches of 
a given symmetry, and we include coupling to the elas- 
tic strains, with e t and r\i the strain component and its 
coupling to r,;. If the phonon Hamiltonian is written 



(12) 



with M an appropriate ionic mass, then the linear cou- 
pling to Qi can be replaced by an effective electron- 
electron coupling by the displaced oscillator transforma- 
tion: 



vT(-q)n{q) 



(13) 



This transformation approximately decouples the 
phonons, leaving an interaction 

Hi n t = -^2vi«in(i) - ]^2Jt{q)n{q)T t {-q), (14) 

i,l «><? 

with Ji(q) = j: m Kr(q), 



(15) 



G 



This decoupling is only approximate, since the Q's de- 
pend on the r's, so that they do not obey canonical com- 
mutation relations. It will, however, be adequate for the 
present purposes. At mean field level, 



iji < a >= m < ri >, 



(16) 



with /ij— j= iV?7?/2c° and cf the bare elastic constant. 
LettingN 



A; 



£(*r(o) -i^jq"®), (i 



7) 



and assuming only one Aj is non- vanishing, the gap equa- 
tions are 



Xi < n >, 



(18) 



< Ti >=-VA i «I«, (19) 
f ■ E + (k)-E^k) 

with E± the energy eigenvalues and f(E) the Fermi func- 
tion. If Ai ^ A2, only the mode with larger Aj has nonzero 
expectation value, although both modes can soften above 
the transition temperature. However, when 



Ai — A2 



(20) 



then Ai = Acos#o, A 2 = Asin# : with 6 arbitrary. 
Equation |2TJ is the generalized version of condition (ii). 



IV. APPLICATION TO THE CUPRATES 

A. Bond Stretching Modes in LSCO 
Three-Band Model 

In Section V, we will calculate the mean-field proper- 
ties of the various CDW's, combining electron-electron 
and electron-phonon interaction and including strong 
coupling corrections and competition with the Mott gap. 
Here we develop a simple model for estimating the 
electron-phonon interaction. We introduce a three-band 
model, and apply it to a number of phonon modes which 
have been proposed to play a role in the cuprates. In 
the following subsection, we will reduce the model to a 
one-band version. 

We concentrate predominantly on CuC>2 plane phonons 
associated with in-plane vibrations, and particularly on 
Cu-0 bond stretching phonons. Mott's original picture 
of a metal-insulator transition involved a loss of cova- 
lency as a lattice is gradually expanded. Hence, near a 
Mott transition it would be expected that bond-length 
changing phonons might play an important role__modu- 
lating a crossover from covalent to ionic behaviorEil. Fig- 
ure H shows the phonon dispersion in La2Cu04 alcpg 
the (C,CjO) (£) direction, calculated by Wang, et al.E3. 



Shown in bold are the branches which are predominantly 
associated with in-plane Cu-0 vibrations. (Wang, et 
al. provide the wave functions only along the high sym- 
metry axes, r and X; the remainder of the curves are 
an interpolation based on symmetry. Note the anti- 
crossing behavior evident on several branches.) The pat- 
tern of distortion for a number of_modes of interest 
also shown, including the breathiru£3, half-breathitigHj'EJ 
(at (|L_0)), dimerfSl, quadrupolarli3, ferroelectric^, and 
shearEj modes. This by no means exhausts all relevant 
modes. Among modes which are not modeled are the 
LTT tilt mode and various c-axis modesE-X As discussed 
above, we treat the LTT mode as a secondary modifica- 
tion of the shear mode instability. 




FIG. 5. Phonon modes of the Cu02 planes, includ- 
ing distortion patterns for several Cu-O bond stretching 
modes. Clockwise from top left: half-breathing, breathing, 
quadrupole, dimer I, dimer II, shear, and ferroelectric modes. 
Filled squares (circles) = Ei (S3) branches predominantly 
representing the Cu02 planes. 

In order to estimate how strongly these modes couple 
to the electrons, we analyze a simple three-band disper- 
sion for the Cu02 plane, with parameters to (t[) for Cu-0 
(O-O) hopping, and A (Ai) for Cu (O) site energy. The 
electron-phonon coupling is assumed to arise via modu- 
lating these parameters. Thus, a change 5u in the Cu-0 
bondlength a produces a change 



St = —[3to5u/a 



in to and a change 



So = ctdAoSu/a 



(21) 



(22) 



in A (or Ai). We estimate ~ 3.5a£30 (a* = 1) and 
ad from a Madelung contribution as ay = qe 2 / eoaAo ~ 
0.55, for Ao ~ AeV, effective charge q = 2 (0 2 ~ or 
Cu 2+ ) and dielectric constant eo ~ 5. In the following, 
we will treat the a's as dimensionless parameters of or- 
der unity in order to scale the electron-phonon coupling. 
The net change in A on any site is then (n — m)So where 
n (m) is the number of near neighbors that move closer 
to (further away from) the given site. For a frozen lattice 
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distortion of amplitude ±5u, we calculate how the net 
electronic energy changes. Model parameters include the 
position of the Fermi level, the value of t[, and the ratio 
[3 /aa- This last parameter can be thought of as a mea- 
sure of covalency: increasing A tends to localize holes 
on Cu, making the material more ionic, while increas- 
ing to delocalizes holes, increasing covalency. We restrict 
ourselves here to the following parameters: to — 1.3eV, 
A = 4eF, a = 1.9 A, Fermi level at VHS (in the ab- 
sence of phonon distortion) at a doping x = 0.25 (to 
achieve this, a large value of tfy/to = —1.8 is assumed); 
and aa — at = 1. The model neglects possible cou- 
pling to the Cu d z 2 orbital (important for the conven- 
tional —LT distortions of the quadrupole and breathing 
modeal3) and anomalies associated with the near insta- 
bility of 2 ~ (Section V.C)£3 The two forms of electron- 
phonon coupling are similar to those assumed by Yone- 
mitsu, et alE3; while the numerical estimates are not 
identical, we find similar estimates for the coupling to 
the breathing mode. Their electron-phonon coupling pa- 
rameters are, in our notation A Q = /3 2 t /Ai, r ~ 0.264, 
X = a 2 d Aj/t A br ~ 0.204a;* = 0.051, for a d = CUa, and 
A br = 60.4ey from Table III. Yonemitsu, et al£3 esti- 
mated A Q = 0.28 and had no estimate for \p. 

For the shear modes and the modes at S = (n,ir), 
the resulting electronic hamiltonian can be written in the 
form H = 



(A + 5 1 


-tie x 


-t2e y 





t 3 e* x 


Ue* y \ 


-t* ie * 


62 


-2t[c' 







-2t[c' + 




-2t[c'_ 


S3 


t(j€y 


-2t' lC ' + 








t* p* 




A -5! 


~he x 


-t%e v 


t 3 6 x 





-2t' lC ' + 




Si 


-2t' x d 


\ ^4^y 


-2t' lC ' + 





—t*p* 


-2t\d_ 


-h I 



(23) 

while for the half breathing mode at X = (it, 0), and the 
ferroelectric mode, H = 



with other modes, the value of S3 will be set to zero. The 
remaining elements are listed in Table I. 

Table I: Phonon Coupling Parameters 



Mode 


Til 


"2 


"3 


mi 


m 2 


m 3 


777,4 


1. Breathing 


4 








1 


1 


1 


1 


2. Half-Breathing 


2 





2' 


1 




1 




3. Quadrupolar 











1 


-1 


1 


-1 


4. Dimcr I 





-2 





1 





-1 





5. Dimcr II 





-2 


2 


1 


-1 


-1 


1 


6. Dimcr III 





-2 


2 


1 


1 


-1 


-1 


7. Shear 





-2 


2 


1 


-1 


1 


-1 


8. Ferroelectric 











1 




-1 





Table II: Energy Parameters 



Mode 


7 


Xi 


x 2 


1. Breathing 





2X\o 


X- 


2. Half-Breathing 





X\o 




3. Quadrupolar 










4. Dimer I 


1 





X + + Ait St s x 


5. Dimer II 








X + + 4ito5to(s x - s y ) 


6. Dimer III 








X + + Ait St (s x + s y ) 


7. Shear 








2 ^ ~fj j I ~fj J 


8. Ferroelectric 








2(tg - Stye* 


9. Flux Dimcr I 


1 





X- + 4dtoStoc x 


10. Flux Dimcr II 








X_ + Mt 5t {c x - c y ) 


11. Flux Dimcr III 








X- +4it St a (c x + c y ) 


12. Flux Shear 








X- - 4t St (s x - Sy) 



The above calculations refer to static (frozen) phonon 
modes, and hence do not give direct information on dJT 
modes. In a Hartree-Fock calculation, dynamic modes 
can be modelled by imaginary order parameters, and in 
Section V, it is shown that the mean-field decomposition 
of nearest-neighbor Coulomb repulsion produces a direct 
coupling to the OAF (flux phase) mode. In this spirit a 
similar mean-field phonon flux phase can be constructed 



(A + 5 l 

-tie* 
2it s' y 


t$e x 





Here 



-t\e x 


2it 1 s y e* 
t 5 e* 




^UoS'y 

-2itiSy€ x 

S 3 


''Zit'is'ye* 






t 3 e* x 


t<se x 








2it^Sye x 


A — 61 




-t 7 e* 





2it S' y 


2it 1 s y e x 



(24) 



= exp(ikia/2), i = x,y, c'± = cos (k x ± k v )a/2 
and s' y — smk y a/2. For the modes at S (the shear 
mode), Si+2 = Si (—Si) i = 2,3. Writing Si = riiSo and 
U = to + rriiSto, then for oxygen vibrations (the breath- 
ing, half-breathing, ferroelectric, and quadrupole modes) 
= +m„ while for modes involving copper motions 
(the shear and dimers), m i+ 4 = — m,*, i= 1,4. For com- 
pleteness, we include an extra parameter ^3 in Eq. |24| 
which is allowed by symmetry, due to second-neighbor 
relative motion; however, to make a uniform comparison 



2it x s y e* 



-2lt Q Sy 

from the above by (a) replacing Sto by iSto and (b) mak- 
ing a similar replacement for the Si, which involves mov- 
ing the term to the off-diagonal (Hij+3) position. 

The lattice distortions are found by minimizing the 
sum of the hamiltonian Eq. ^ and a phonon term, 
Eq. |l2], which may be rewritten as H p h = A(5u/a) 2 , 
where A = r]Muj 2 a 2 /2, and A and u> are listed in Table 
III. Strictly speaking, the phonon frequencies should be 
bare values, neglecting electron-phopon coupling, but foe 
present purposes we use calculated^ or experimental 
values for the undoped insulator La2Cu04. [For the shear 
mode, A = (Cn-C 12 )a 2 c, where {C u -C 12 )/2 = 99GPa 
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[ |36| and c is the c-axis lattice constant.] The parameter 
r\ is one if the distortion is along the X or Y axis only, 
two if it is along both, and M is the mass of the moving 
ion, oxygen or copper (or the appropriate average for the 
ferroelectric mode). These equations must in general be 
solved numerically, but some insight can be gotten for 
the special case n 2 = n$ = t\ = 0, for which the four 
non-zero eigenvalues are given by 



E i± = 



Ap ± y/Al + AW 2 



(25) 



withi = ±1, Wg = W?±y/X'f + \X 2 \ 2 , Wf = 4(t§+(l- 
•y/2)6t 2 ) (except for the half breathing and FE modes, for 
which Wl = 2(tg[l + 2s'*} + StD) and the X'a and 7's 
are listed in Table II, with Si = sin(kia) (and recall Cj = 
cos{k l a)) X± = 2{t 2 ± 8t1)(c x + c y ) =p ^l8t\c y and X w = 
2E8\ + itoSto (note that the solution is only implicit if 
5i ^ 0). The generic form of X 2 is X 2 = t 1 t 5 el+tltfe* x 2 + 
ht&e 2 + title* 2 . Note that E l+ is the antibonding band 
which is approximately half filled. 

Figure || shows the calculated electronic energy shifts 
up to large distortions Su/a = 0.1 (local distortions of-up 
to about 5% are reported in EXAFS measurementscj) . 
A number of interesting features can be noted. (1) All 
the distortions lead to highly nonlinear shifts of the elec- 
tronic energy. (2) While the breathing mode has [the 
strongest coupling (as originally predicted by WeberEZI), 
the half-breathing mode also has a large coupling, since 
Si produces a large splitting at the VHS, Fig. ||. (3) The 
energy shifts correspond to significant phonon softening, 
but when the quadratic phonon energy is included, it 
is found that only the breathing and (marginally) half 
breathing mode become unstable - the additional soft- 
ening due to electron-electron coupling is discussed in 
SectioiuV. Instabilities of the dimer modes were found 
earlierEj, but only in strong nesting conditions - half fill- 
ing with t[ = 0. (4) A problem for a pure phonon- flux 
phase should be noted: the symmetric dimer III mode has 
a larger energy lowering than the asymmetric (c x — c y ) 
dimer II mode. 




FIG. 6. Energy of various phonon modes of CUO2 planes, 
assuming = 1 and (3 = 3.5, for static (solid lines) or flux 
phases (dashed lines). The phonon modes are numbered ac- 
cording to Table II. 




FIG. 7. Doping dependance of the frozen phonon en- 
ergy lowerings, associated with the breathing mode (cir- 
cles), half-breathing (n/a,0) mode (squares), shear mode (di- 
amonds), and dimer III mode (triangles), all normalized to 
their maximum amplitude. The half-breathing and shear 
modes have strongest coupling near the VHS. Su/a = 0.04 
for the breathing mode, 0.08 for the others. 




FIG. 8. Three-band dispersion of the CUO2 planes, in 
the presence of a frozen breathing mode (dashed lines) or 
half-breathing (jr/a, 0) mode (solid lines) modulation, assum- 
ing P/a d = 3.5, and Su/a = 0.02 (breathing) or 0.04 (half 
breathing). The arrow indicates the breathing mode hole 
pocket near (n/2a, it /2a). Brillouin zone special points are 
r = (0, 0), X = (tt/o, 0), S = (tt/o, ir/a). 

The doping dependence of the energy lowering is shown 
in Figure pi— ■ The results are in good agreement with 
experimenfapJ on both LSCO and YBCO and with LDA 
calculations^ which find that the half-breathing mode 
softens more than the (it, ir) breathing mode. This re- 
sult is surprising since the breathing mode has stronger 
coupling. However, this mode has an unusual frustration 
effect: it opens a complete gap at the Fermi level, Fig. ||, 
so the Fermi level can only fall in the gap if the band is 
filled, x = 0. To accomodate additional holes, the Fermi 
level must shift below the gap, reducing the electronic en- 
ergy lowering. Hence optimal doping for this mode is at 
half filling, and while there is a large electronic softening 
of the breathing mode at the VHS, there is an even larger 
softening near half filling, so a doping dependence would 
show a hardening of the mode. Indeed, such a hardening 
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is observed in the nickelatesii. This result persists in the 
presence of strong correlation effects, Section V. For the 
half-breathing mode correlation effects should enhance 
the already large softening of the mode with doping, in 
good agreement with experiment. 



\ 

\ 

\ 

\ 


y? / — «A 




< / i \ 








^\ / \ 


/ 7 






v 






\ 





X,Y 



FIG. 9. Three-band dispersion of the CuC>2 planes, in the 
presence of a frozen dimer II mode (long-dashed lines) or shear 
mode (solid lines) modulation, assuming 13/ ad = 3.5, and 
du/a = 0.04. 



Figures |8|-|10| compare the dispersions of several of the 
modes, near the upper (Cu-like) band of the three-band 
model. For all cases, the magnitude of electronic energy 
lowering can be directly correlated with the size of the 
gap at the VHS (X-point). Note that the gap is largest 
for the breathing mode, but the Fermi level lies below 
the gap if x ^ 0. For the modes at T, there is no unit cell 
doubling, but the dispersion is different along X and Y, 
leading to an effective splitting of the VHS's, Figs. |9-10 

For the half-breathing mode we find that the softening 
is via the on-site Cu energy, A. Making A inequivalcnt 
on alternative rows of Cu's along the X-axis (parallel to 
the oxygen distortion) leads to a large splitting of the 
VHS's. Figure || shows that the gap is smaller than for 
the full breathing mode, but is clearly centered at the 
VHS. 




X,Y 



FIG. 10. Three-band dispersion of the Cu02 planes, in the 
presence of a frozen ferroelectric (solid lines) or quadrupole 
mode (dashed lines) modulation, assuming /3/ay = 3.5, and 
Su/a = 0.08. 



B. Bond Stretching Modes in LSCO 
One-Band Model 

To include electron-electron coupling, it is convenient 
to first reduce the above results to a one-band model. 
The parameters of the one-band model can readily be re- 
lated to those of the three band model by expanding the 
approximate eigenvalues of Eq. as Ej + = Ao+Wf/Ao. 
In this way, it is found that the effective Cu-Cu hopping 
parameter is t = t\j Ao [with the parameters assumed 
above, this would give t = 0A2eV, somewhat larger 
than the value 0.326eV estimated from photoemission] , 
with phonon distortion parameters St/t = 26to/to and 
Si = Sq + St. Linearizing the results in Sto, it is read- 
ily seen that the flux dimer II mode has the symmetry 
of the conventional flux phase, and it is found that the 
quadrupole and ferroelectric modes do not couple at the 
one band level, while the breathing mode couples only 
via S\. The electronic dispersion can be solved: 



E = -4t'c x c y ± J{m 1 ~5 1 f + \X X + X 2 



(26) 



with mi and X 2 listed in Table III, and X\ = — 2t(c x 
c y ). The result for the half breathing mode is distinct: 



E = -2tc y ± 



(2£i 



2Stc y ) 2 + A{t + 2t'c y ) 2 cl. (27) 



The resulting dispersions are in excellent agreement with 
the corresponding antibonding band dispersions of the 
three band model. 

For the half breathing mode the term in St follows from 
the ^3 term in Eq. ^J: for a phonon distortion along X , 
the hopping along Y is modified along alternate rows, 
at t ~ £q/(Ao ± ^3). Such, a term has been previously 
considered by Shen, et al£3. In the present comparisons 
of different modes, the role of such a term has been ne- 
glected (by setting ^3 = 0). 

Table III: One Band Model Parameters 



Mode 


mi 


x 2 


w(mcV) 


A(eV) 


1. Breathing 


4 





66 


60.4 


2. Half-Breathing 


2 




83 


47.8 


3. Quadrupolar 








56 


43.5 


4. Dimer I 





—2iSts x 


38 


39.7 


5. Dimer II 





—2iSt[s x + s y ) 


38 


79.5 


6. Dimer III 





—2i8t(s x — s y ) 


38 


79.5 


7. Shear 





-2St(c x - Cy) 




59 


8. Ferroelectric 








77 


32.8 


9. Flux Dimer I 





— 2iStc x 


38 


39.7 


10. Flux Dimer II 





~2iSt(c x - c y ) 


38 


79.5 


11. Flux Dimer III 





—2iSt(c x + Cy) 


38 


79.5 


12. Flux Shear 





+2iSt(s x — s y ) 




59 



V. RESULTS IN STRONG COUPLING LIMIT 
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A. Electron-Electron Coupling 

In one-dimensional metals, the CDW instability in- 
volves both electron-electron and electron-phonon cou- 
pling. Here we show that the same is true in the present 
two-dimensional problem: electron-electron interaction 
actually dominates the electron-phonon coupling in the 
cuprates. A nearest neighbor Coulomb interaction V con- 
tributes to all threes-components of pseudospin, Eqns. ||- 
pj; it is estimated^ that V ~ 0.2 - 0.3eV. 

The electronic part of the Hamiltonian can be written 
H = H e j + Hi where the electronic kinetic energy is 



Ei E 
E Ei 



(28) 



with E = —2t(c x + c y ), Ei = —kt'c x c v , and t (f) is the 
(next) nearest neighbor hopping parameter. At mean 
field, the effects of both V and the on-site Coulomb in- 
teraction U can be included, leading to the Hamiltonian 
matrix 



Note that while the interaction Hi resembles H v u, of 
Eq. [27], there are a number of differences. First, while 
the term involving T y is generated dynamically in the 
molecular problem, here a term in R y appears explicitly. 
On a single plaquette, the three solutions are degener- 
ate (the eigenvalues depend only on R, and not on the 
individual i?i's), but intercell hopping and fc-dependent 
coupling lead to slight dependence of the dispersion on 
the individual Rio 's. (Note that R x simply renormalizes t 
to have unequal values in the x and y directions—The ef- 
fects of this term have been recently analyzedc3o.) The 
solutions remain degenerate at the VHS, for equal R^s. 

Introducing the polarization 



k.t 



E A 



E_ 



the gap equations can be written 



E [^,^ = 1, 

- J k 



(31) 



(32) 



Rza ^ 
Eq + R x 



Ei 

— iR y 



Eq + Rx 
— Rz.a 



- iR y 
■ Ei 



with eigenvalues 



E 



±,<7 



= E X ± JR 2 



El 



2R x Ei 



0, 



(29) 



(30) 



where R 2 = R 2 X + R 2 y + R 2 Z(J , Fig. y. In Eqs. 
we take R x = Rxol, Ry = Ryol, and R ZjCr = R z0 + 
aUm z , where the i?io's are constants, a = ±1 is the 
spin, and m z is the average magnetization per site in 
the antiferromagnetic phase induced by U. Parameter 
values are given in Section MA. The conventionalEI flux 
phase dispersion is recovered from Eq.|3(] for the special 
choice of parameters: t' = R x o = R z , a = and R y o = 2t 
(compare Eq. H). 
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for i including all the non-vanishing gap contribu- 
tions. The most general case has all four gap pa- 
rameters, R x q, R y o, R z o, m z , nonvanishing, with corre- 
sponding Ai = Ai(7 + Eq/R xQ )j, A 2 = A27 2 , A 3 = 
h{Rz,a/Rzo), and A 4 = A 4 (l + aR-<JUm z ). For the 
extended Hubbard model, the A's areE3: Ai = A2 = 8V, 
A3 = 16V — 2U, and A4 = 2U. (In the areas of overlap, 
these gap equations agree with those found by NayakEl 
However, a strong coupling calculation^ has suggested 
that only an attractive near-neighbor coupling would sta- 
bilize the flux phase.) 

In the absence of on-site U, the near-neighbor inter- 
action V favors the CDW, R z . The energy lowering is 
exactly twice that of the other components, the factor 
of two arising from spin (since the x and y terms in- 
volve hopping, only terms with the same spin on both 
neighbors contribute; the z term has no such limitation). 
This is counteracted by strong coupling effects, since the 
CDW involves an imbalance of Cu site occupancies. Ap- 
proximating V = t, U = 6i, 8V > 16V — 2U, favoring 
the flux and shear phases. Because the couplings are not 
symmetrical (Eq. |3^) , these modes have slightly different 
gaps, R x o = 340meV, R y o — 3207716^, R z o = 8AmeV. 
However, 2U > 8V, so in a purely electronic model the 
antiferromagnetic phase dominates, both at half filling 
and in the doped case. In the following subsections, it 
will be shown that electron-phonon coupling can reverse 
this situation in the doped cuprates, and also brings the 
CDW energy closer to that of the other phases. 



FIG. 11. Energy disp ersion in the presence of pseu- 
dospin ordering, Eq. |(| for t = 0.326eV, t' = -0.276£, 
and R x o = R y o = RzO = O.leV (short dashed line); or 
Rio = 0.1732eV and the other two R's = 0, for i = x 
(dot-dashed line), y (solid line), or z (long dashed line); or 
all Rt's = (dotted line). 



B. Electron-Phonon Coupling 

Combining the above results with electron-phonon 
coupling, we limit discussion to the three modes of 
primary interest. In the one-band model, the linear 
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electron-phonon coupling has exactly the same form as 
the V electron-electron term in Eq. 29. Thus, the free 



energy F is a function only of the sums of the electron- 
electron and electron-phonon gaps, leading to a simplifi- 
cation of the gap equations. For example, for the CDW 
F = F(R z0 + 46 1 ), so 



^ dF_ ^ OF 
d5i dR z0 



(33) 



Since the averages are just those evaluated in the gap 
equations, they can be replaced by the derivatives of the 
quadratic terms, yielding: 



Rzo 



A(16V-2U) 



0.22 



(34) 



for ad — 0.5 and the other parameters as given above. 
For the shear mode, the equivalent result is 



St 
Rxo 



pH 2 
AV 



0.068, 



(35) 



with at = 1; an identical result holds for the flux phase, 
with appropriate A, yielding 5t/R y o ~ 0.05. These are 
the analogs to Eq. |l3|. They are weak coupling results 
and should be modified when the gaps become compa- 
rable to the phonon frequencies. Thus, (1) there is sub- 
stantial phonon coupling, but the transitions appear to 
be predominantly electronically driven; (2) the largest 
correction is for the CDW mode, although the correction 
is not large enough to make the CDW as unstable as the 
other two modes (if a^ were twice as large as estimated, 
all three modes would be approximately degenerate) . Fi- 
nally, (3) near half filling the corrections are not large 
enough to make any of the structural anomalies compet- 
itive with the AFM instability-JThis result is consistent 
with the finding of Hsu, et al.E3, that the flux phase is 
unstable against magnetic order near half filling. It re- 
mains possible that additional couplings could tip the 
balance toward structural distortions, particularly in the 
doped materials. Certainly, there .are strong indications 
of deviations from Migdal theoryEj, and earlier calcu- 
lations of the-JICO and LTT instabilities found strong 
nonlinearitiesoo. In the following subsection we explore 
another possibility, a particular form of correlated hop- 
ping. 



C. Unconventional Coupling Associated with O 2 

The O 2- ion is known to be inherently unstable, be- 
ing stabilized in a solid by the Madelung potential of 
surrounding ions. This near instability has been sue-, 
gested to be a driving force in ferroelectric transitionsE.il 
in perovskites and in cuprate superconductivityo. One 
of the principal manifestations of this near instability is 
the large change in ionic radius on doping 2 ~ to 0~.E£l 



While this anomaly should affect the electronic proper- 
ties in a number of ways, we will here explore only one 
aspect, a correlated hopping. 

When a hole is localized on a single oxygen, the shrink- 
age of its radius allows the adjacent coppers to approach 
much more closely. This would enhance the correspond- 
ing hopping probability, t, except that the presence of 
the hole inhibits other holes from hopping to the same 
site. If, however, the first hole hops away, it will take 
the local lattice some time to relax back, and in that in- 
terval there will be an enhanced probability of another 
hole hopping onto the oxygen. In a one band model, this 
would correspond to correlated hopping between adja- 
cent Zhang- Rice singlets, or effectively, between adjacent 
coppers. This hopping adds a term 



H c . h . = -t cl c h c j° c }<7' 

<i,j> ,<i' ,j> 



(36) 



to the Jiamiltonian. Such a term has been condidered Jay 
NayakEIll; for related models, see Hirsch and Marsiglioo. 
This term contributes 96t c i to all three coupling parame- 
ters Aj, i = 1, 3. From the present considerations, a value 
t c i — —ZxPiocPtSuo/a ~ xt can be estimated, where 
Suq ~ — 0.3A is the shrinkage in radiusL3 on going from 
2 ~ to 0~ . The factor x arises because holes only move 
preferentially onto the oxygens for doping beyond half 
filling, x > 0, while /3/ oc < 1 is a parameter introduced 
to describe the degree of localization of the hole on a 
single oxygen - predominantly due to polaronic effects. 
This factor is related to the delicate issue of the crossover 
from ionic to covalent behavior: f3i oc = 1 corresponds to 
the ionic limit. The t c \ correction is so large that even a 
significantly covalent correction /3; oc < 1 would allow the 
structural distortions to overcome the AFM state - but 
only for x > 0. 



D. Numerical Results 

Here we describe our numerical results. The gap equa- 
tions follow from minimizing the mean-field free energy 



F = Yl EifiEj-TS + Nfo 



(37) 



where S is the entropy, N the total number of electrons, 
and 

= Rio + Ri + + y, A Sut a + 2 + g 



W 8V-U 



i=l,3 



(38) 



The Ei are solutions of Eq. |30j with phonon coupling in- 
cluded by the substitutions: R z q — > R z q + AS±, R^ — > 
Rio + 2<5£, i = x,y. In Figure [l2] we gradually turn on 
the phonon coupling, by varying the a factors from zero 
(no phonon coupling) to 1, covering the expected range 
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of coupling. Here we study the coupling near the VHS, 
x = 0.25. What is plotted is the self-consistent con- 
densation energy. For weak phonon coupling, the CDW 
has the weakest coupling, due to the on-site Coulomb re- 
pulsion. However, it has the strongest intrinsic coupling 
to phonons, so as ad is increased, it grows fastest, and 
can actually cross the other modes. The dotted line is 
the energy of the simple doped antiferromagnet (also at 
x = 0.25). It can be seen that for reasonable phonon cou- 
pling, the ground state is paramagnetic. The dashed lines 
show that even moderate correlated hopping (/3; oc = 0.1) 
significantly enhances the binding energies, and shifts the 
degeneracy points of the CDW with the other modes to 
lower values of a. 
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FIG. 12. Binding energy of various modes - site CDW (cir- 
cles), bond CDW (squares) and flux phase (diamonds), for 
flioc = (solid lines) or 0.1 (dashed lines); horizontal dot- 
ted line = energy of doped antiferromagnetic phase. For all 
curves x = 0.25. 
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FIG. 13. Doping dependence of binding energies of phonon 
modes: site CDW (circles), bond CDW (squares), and 
flux phase (diamonds), assuming all mode a's = 0.6, and 
Pioc = 0.1(a;/0.25). Also shown is the doped antifsHjomag- 
net (triangles) and ferromagnet (inverted triangles) p3]. The 
dotted line is the tangent construction for phase separation. 
For convenience in viewing, an energy shift 5E — 2.1(1 — x)to 
has been added to all curves. 



In Fig. |13| the doping dependence of the binding ener- 
gies is illustrated for a representative a, = 0.6 (i = d, t), 
compared to-, that of the doped antiferromagnet and 
ferromagnetEZl. All three CDW modes have a parabolic 
binding energy vs doping: for the bond-CDW and flux 
modes, the energy minima are near the VHS doping, 
while the site-CDW minimum is shifted to lower doping 
by the frustration effect discussed above (the minimum 
is no longer at x = 0, due-iq the strong correlation ef- 
fect). A by now familiarliZlo Van Hove induced phase 
separation could arise between the antiferromagnetic in- 
sulator at half filling and any of the CDW modes. For 
the present parameter values, the lowest energy state in- 
volves the bond CDW at x ~ 0.27. Note in particular 
that the boad CDW is more stable than the ferromag- 
netic phaseH3. This is consistent with the numerical re- 
sults of Yonemitsu, et al.Ej, who found a crossover from 
magnetic to dielectric polarons in the low doping regime. 
Remarkably, all three CDW curves approximately con- 
verge (weak JT degeneracy) at a doping, x ~ 0.21, close 
to the VHS. However, this crossover point is not directly 
observable, due to the phase separation. At this time we 
have not made a detailed study of the parameter depen- 
dence of these results, but the bond-CDW appears to be 
most stable. In particular, (1) ad should probably be less 
than at; and (2) our estimate for V may be a little large, 
but t c \ is probably too small, so V can be reduced and 
t c i increased to keep a fixed coupling to the shear mode. 
Both of these corrections would act to further weaken the 
site-CDW. The results are quite intriguing. For a reason- 
able estimate of the parameters there is phase separation 
between the antiferromagnetic insulator and a param- 
agnetic state associated with LTT-like distortions, while 
the CDW and flux phase states are close in energy. 

These results are in good agreement with the numerical 
calculations of Yonemitsu, et al.E3. In their work V was 
not included, so their phonon anomalies are of the breath- 
ing mode type, but they see clear evidence for nanoscale 
phase separation between antiferromagnetic background 
and dielectric polarons. 

The present approach of turning on the electron- 
phonon interaction last (Fig. |l2|) makes the resulting 
near-degeneracy seem rather accidental. It is more nat- 
ural to reexpress the result: for both electron-electron 
and electron-phonon coupling, the flux and bond-CDW 
phases are nearly degenerate, while the site-CDW is more 
strongly coupled. However, the on-site Coulomb re- 
pulsion opposes site-CDW formation, making the three 
modes nearly degenerate in energy. 

VI. DISCUSSION: FUTURE DIRECTIONS 

The present paper establishes the framework for ana- 
lyzing CDW's in the cuprates, and Berryonic matter in a 
number of related materials. Clearly much work remains 
to be done. This includes: 
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(1) From the associated susceptibilities the modifi- 
cations to the electronic disBpsion can be calculated 
and compared with the kinl£3'E2l seen in photoemissioa. 
While-phononic contributions appear to be important^, 
Vallate has noted that the linear frequency dependence 
of the imaginary self energy extends to much too high fre- 
quencies to be solely due to strong electron-phonon cou- 
pling. This is consistent with the present results, that the 
CDW has a strong electron-electron component. Hence, 
both components will be important in interpreting the 
photoemission dispersion, with proxiewty to a VHS pro- 
viding a marginal Fermi liquid-like£23 background and 
phonon coupling the kink in the dispersion. 

(2) The manner in which the y-component couples to 
the Hamiltonian is markedly different in the molecular 
and lattice versions of the theory: in the molecular ver- 
sion, the direct electron-phonon coupling is to the x and 
z components {B\ and E>2 modes), with the dJT term 
generated dynamically - from the phonon kinetic energy 
term. In the lattice, the flux phase term arises directly 
from electron-electron coupling via V, and corresponding 
terms are allowed in the Hartree-Fock expansion of the 
electron-phonon coupling. It remains to be seen what 
role phonon kinetic energy plays in the lattice problem. 
[One effect will presumably be the kink in the electronic 
dispersion when E — Tiujph-} 

(3) Fluctuations of the CDW's (or of the stripes) 
should apuaar as new low frequency (or pinned) pha- 
son modesL3. It is possible that such-modes have been 
observed in microwave measurements EUd 

(4) It remains to work out the detailed doping depen- 
dence, in terms of stripes. 

(5) It will be important to provide detailed calculations 
showing how the shear mode couples to the various bond 
bending modes in the different cuprates. A related issue: 
are there tilt distortions in the flux phase? (if so, then 
experimental evidence for local LTT order could really 
be associated with the flux phase.) 

(6) One puzzling feature is the role of the LTO phase. 
Near optimal doping it would seem to be predominantly 
associated with local LTT order, but it seems to be a 
uniform phase near half filling. In principle, its presence 
could be accidental, particularly since it is only present in 
LSCO, but it does provide an easy axis for orienting the 
spins, and hence may be involved in the stripe crossover 
from vertical to diagonaL (Also, the pseudogap seems to 
follow T pg (x) = 2T LTO B) 

(7) A very similar model should apply to the nickelates, 
and extensions can be made to other forms of Berryonic 
matter. 

(8) There remains the problem of the competition of 
density wave order with superconductivity, and the pos- 
sible roles of magnetic fluctuations. 



VII. CONCLUSIONS 

This results of this paper have bearings on four sep- 
arate issues. First, on a purely formal level, the paper 
presents an alternative approach (phonon flux phases) 
to the dJT effect in solids, revealing previously unsus- 
pected analogies between flux phases and Berryonic mat- 
ter. Secondly, The approach is applied to the cuprates, 
which are an a priori unlikely candidate for JT effects. 
Nevertheless, a 'hidden' JT degeneracy is found, and a 
numerical estimate suggests that the cuprates are close 
to the weak degeneracy point of a square molecule, which 
should enhance the possibility of a dynamic flux phase. 
These calculations can be considered a generalization of 
the results of Yonemitsu, et al.E3, by (a) reducing the 
problem to a one band model, (b) including a nearest 
neighbor Coulomb repulsion which enhances the scope 
of CDW- like instabilities, and (c) providing a mean- field 
underpinning for the numerical calculations of local elec- 
tronic phase separation. Third, the paper explores the 
relative strength of different CDW-like distortions in the 
cuprates, and suggests an explanation for strong coupling 
to the LTT and half-breathing modes. The role of these 
distortions in stabilizing charged stripes is discussed in 
an Appendix. Finally, the results offer strong support 
for the general picture of stripe phases as stabilized by 
VHS-induced ordering. A similar approach should be 
applicable to other systems, and in particular, to nicke- 
lates where the coupling to the CDW/breathing mode is 
known to be stronger. 

For the cuprates, we find that the structural distortions 
are close to the dJT degeneracy, but that the shear and 
flux phases are more unstable than the CDW. Hence, 
one would expect significant softening of the breathing 
modes, but instability in either the shear or flux phases. 
These expectations are borne out in experiments on the 
cuprates, with the shear mode coupling to a (local) LTT 
order. The near degeneracy of the shear and flux phases 
is consistent with experimental evidence for both modes. 
Perhaps the most interesting phase is the flux phase, 
which is closely related to a chaotic dJT phase. 

Recently, dJT phases have been proposed in a numbec 
of exotic materials, including cuprateaij, buckyballsLHl 
and manganatesES. To describe the dJT phase of buckjj.- 
balls, the concept of Berryonic matter was introduced!! 
In this model, no attempt was made to accurately model 
Buckyball solids. Instead, a lattice of dJT molecules was 
assumed, basing their properties on the known anoma- 
lies associated with triatomic molecules. Since the square 
molecule, Fig. 0, has the same dJT anomaly as the trian- 
gular model when Ej T = Ej T , and this degeneracy per- 
sists on a square lattice (Section IIIC), the present model 
should be an excellent starting point for studies of Berry- 
onic matter. Moreover, the proximity of the cuprates to 
the Berryonic limit should stimulate interest in this un- 
usual state of matter. 
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APPENDIX A: PHONONS ON STRIPES 

In the Van Hove model, charged stripes (or domains) 
are stabilized at a particular density by.-a VHS driven 
ordering instability, most likely a CDWliU as analyzed 
here. Originally, this distinguished the VHS stripe model 
from magnetic models, in which the charged stripes were 
structureless domain walls between anti£ecromagnetic do- 
mains. However, in 1995, it was foundE23 that the hole 
density on the charged stripes in LSCO:Nd is consider- 
ably less than one per copper, and the magnetic mod- 
els weij£_mo dined to include CDW order on the charged 
stripesEHj. The form of CDW order proposed gener- 
ally differs from that discussed here (being more one- 
dimensional), and the models still do not explain the 
experimental observation that charge ordering arises at 
higher temperatures that spin ordering. Here we dis- 
cuss experimental evidence linking the phonon anomalies 
specifically to charged stripes. It should be noted that 
there must in general be a coupling to all modes: the 
charge modulation will act to scatter all phonons; here 
we are speaking of a stronger coupling, with phonons lo- 
calized on the stripes. 

It is by now clear that the LTT phase plays an impor- 
tant role in stabilizing a nearly static stripe phase£22l in 
La 2 _ a; _j / R y A a ;Cu04 (LACO:R) with R a rare earth (RE), 
typically Nd or Eu, and A = Sr, Ba, or Ca. The reasons 
for this stabilization are less clear. One explanation 
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a simple pinning effect: the stripes 'prefer' to run paral- 
lel to the Cu-O-Cu bonds, and the LTT distortions lock 
them into an ordered array. However, it is not necessarily 
true that LTT order enhances stripe pinning: in Eu sub- 
stituted samples, the Cu spies, appear to be even more 
dynamic in the LTT phasetM Moreover, such a pic- 
ture completely fails to explain the common occurence 
of simultaneous fluctuating stripe order and fluctuating 
LTT order. States with fluctuating local order are rela- 
tively uncommon, and for two dynamic phases to appear 
simultaneously and to interact strongly, but to have com- 
pletely independent origins stretches credulity. It is much 
more likely that the connection is fundamental: that the 
structural and stripe orders are two aspects of a single 
phenomenon. Indeed, this must be the case: neutrons are 
insensitive to individual electrons, so the neutron diffrac- 
tion evidence for long range stripe ordert2HI is really de- 
tecting an accompanying structural order. Perhaps part 
of the problem is a limited vocabulary for dealing with 
fluctuating local structural order: the average macro- 
scopic structure will generally appear to be LTO or LTT 
or some intermediate structure; only if the domains are 
of sufficiently large scale would one expect to see a clear 
mixture of two phases. 

Hence, this Appendix deals specifically with evidence 
that the structural order is also short range, and evolves 
with doping in parallel with the electronic stripes. The 
cumulative evidence supports a picture of coexisting 
structural and electronic stripe structure: whenever 



stripes are present, not only are the electronic proper- 
ties inhomogeneous (with insulating magnetic and con- 
ducting charged stripes), but also the lattice distortions, 
with local LTT distortions on the charged stripes, while 
low-tempcraturc orthorhombic (LTO) distortions are fa- 
vored by the magnetic stripes. This strongly supports a 
model of the stripe phase as a form of nanoscale phase 
separation, with th e . charged-stripe phase stabilized by 
an LTT-like phaseBEl 



1. Coupling of Charged Stripes to local LTT Order 

a. LTT-LTO Competition 

Uniform, macroscopic LTT order is rare in the 
La2- a; Sr a ;Cu04 (LSCO) family of cuprates. There is 
a narrow regime of LTT phase near optimal doping in 
La2_ a; Ba 2 ;Cu04 (LBCO), while long-range LTT order 
can only be induced in LSCO by replacing some La with 
a RE, typically Nd or Eu. However, short range LTT 
order appears to be more pervasive. 

EXAFS measurements by Bianconi and coworkers 
found evidence for anomalcps-long Cu-0 hands in doped 
Ba 2 Sr 2 CaCu 2 8 (BSCCO)H and LSCOB, which were 
interpreted in terms of stripes with aiteipating LTO-LTT 
orders. A number of experimentsc2rca have explored 
the possibilitjiij that the high-temperature tetragonal 
(HTT) and LTO phases can be 'dynamic LIT' phases 
(i.e., with local LTT order). Billinge, et al.c£l analyzed 
neutron diffraction pair distribution functions (PDF's) 
of LBCO, and found that in the LTO phase the local 
order has LTT symmetry, crossing over to average LTO 
structure on a 1(M length scale. Moreover, there is no 
change in the local structure when waxming into the HTT 
phase- EXAFS results are consistento. In LSCO, Bozin, 
et al.ErC-a find a crossover from local LTO order very near 
half filling, to a clear mixing of local LTO and LTT or- 
der in the doped materials. The results are influenced by 
large local distortions near the dopants, particularly in 
LBCO associated with the. large difference in ionic size 
between Ba +2 and La +3 ,cil hu±-.also near the dopant Sr 
sitesEJ. It has been suggestec£23 that the distortions are 
characteristic of the hole-doped Cu0 2 planes, but with 
enhanced distortion near Sr due to partial hole localiza- 
tion. These doped materials display considerable tilt and 
bond-length disorder, interpretable in terms of stripes. 

A crossover from local LTO to local LTT order with in- 
creased doping is an intriguing possibility, particularly in 
light of the crossover in stripe orientation from diagonal 
at low Hopin g! 107 ! which could be 'pinned' by LTO order, 
to horizontal near optimal doping. (It has been suggested 
that the .LIT distortions orient the stripes along the Cu- 
O bondsli23.) The connection between the LTO phase 
and magnetic stripes may be related to strong magnetoe- 
lastic coupling: in the Neel phase, the-spin flop transition 
is strongly suppressed by LTT ordered, while in the diag- 
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onal stripe phase the cluster spin glass transition can be 
seen in anelastic relaxation measurementst£3. The Mott 
insulator S^CuC^C^ (SCOC) is found by x-rays to be 
tetragonal, but below the Neel temperature the infrared 
absorption peak associated with the Cu-0 bending mode 
phonon is found to split into two components^ 



b. LTT and Stripes 

Recent studies have added considerable evidence for 
direct phonon coupling to stripes. First, of course, 
direct neutron diffraction evidence for charge orderE23l 
what is seen by neutrons is not the electronic charges 
themselves, but an accompanying lattice distortion, sug- 
gesting that the distortion is a fundamental part of the 
stripes. These anomalies, and related anomalies of the 
half-breathing modeE3 need not imply—that stripes are 
stabilized by electron-phonon couplingtll: it could sim- 
ply be that there is stronger electron-phonon coupling 
on_hole-doped stripes. On the other hand, Bianconi, et 
al.Ej report EXAFS evidence for a quantum critical point 
(QCP) - a local splitting of the Cu-0 bond length which 
can be tuned to zero by adjusting the chemical micros- 
train. It is found that the highest critical temperatures 
are associated with this optimal degree of 'microstrain'cJ. 
The associated bondlength splittings are much larger 
than can be accounted for by the doping dependence of 
the average lattice constant, and hence point to a struc- 
tural QCP distinct from stripe formation. Finally, there 
is the fmdingC-3 that local LTT-LTO type structural dis- 
order is characteristic of the stripe regime, and the struc- 
tural disorder greatly reduces near x = [ft^^-close to the 
point where the stripe phase terminatestlJ'tiHl. This near 
coincidence of stripe disorder and structural disorder pro- 
vides very strong evidence that the two kinds of stripe are 
associated with two kinds of structural order: that the 
LTO phase is connected with magnetic stripes, at half 
filling, and the LTT phase with charged stripes, near op- 
timal doping. The data cannot be understood in terms 
of preexisting stripes accidentally pinned by a uniform 
LTT order. 

Thermal conductivity measurements suggest a very 
close connection between stripe order and the LTT phase. 
Thermal conductivity k has proven to be an important 
probe of thp-sharge ordering transition in nickelates and 
manganitealid. Typically, k is suppressed above the tran- 
sition due to collisions be 
stripes. Baberski, et al.tli 



ween phonons and fluctuating 
find a very similar suppression 
of k in the cuprates, but always at the LTT transition. 
Figure [l4] compares the thermal conductivity near the 
LTT transition in LSCO:Nd[ [nj] (b) with that of the 
charge ordering transition in a manganiteELj (a), suggest- 
ing an intimate connection between the LTT and charge 
ordering transitions. In the absence of long range LTT 
order, k is suppressed over the full temperature range - 
this is found both for pure LSCO and for the RE substi- 



tuted materials, for x > 0.17. The thermal conductivity 
suppression persists down to the lowest doping studied, 
x = 0.05, and the jump at Tltt actually increases with 
decreasing x. However, suppression is absent in the un- 
doped materials. It would be most interesting to see if 
the suppression terminates near x — 0.05, where there 
is a superconductor-insulator transition and the stripes 
cross over from horizontal to diagonal. 
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The close connection between the LTT phase and 
charged stripes is further illustrated by . g .. repent 
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study of the rich phase diagram 
Lai.6- a; Ndo.4oSr ;I .Cu04, Fig. |l~5|a. This complex behav- 
ior can be simply understood (Fig. |i~5|b) if it is kept in 
mind that the structural phases are macroscopic aver- 
ages, while local structure can be more complicated. We 
assume that the LTO phase is associated both with the 
high-temperature electronically disordered phase, and 
with the magnetically ordered phase at x = 0, while 
the LTT phase is characteristic of hole-doped stripes. 
As temperature is decreased, long-range stripe order 
arises in two steps: first the charges order (filled circles) 
then the spins (dot-dashed line). The other curves in 
the diagram are associated with fluctuating stripe order. 
Thus, the NQR 'knockout' line may indicate the onset of 
stripe fluctuations 



]; atjaearly the same phase bound- 
ary XANESlla and PDFe.! measurements find the onset 
of local Cu-0 bondlength variations suggestive of local 
LTT-type tilts. Near 1/8 doping (vertical line) there 
is a crossover in the stripe fluctuations: at low doping 
there are narrow charged stripes fluctuating in a mag- 
netic (LTO) background; at higher doping narrow mag- 
netic stripes fluctuate in a charged (LTT) background. 
Long range charge order (filled circles) coincides with the 
lower of the LTT(-like) transition or the NQR knockout 
line. 

Magnetic field— studies provide additional evidence. 
The Hall effeciti3 in strongly Nd substituted LSCO 
shows a crossover from one-dimensional (Rh -> as 
T — ► 0) to two-dimensional behavior at x = 0.12; for 
x < 0.12, the Rh anomaly turns on with stripe order 
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- exactly at the LTT-like phase transition. Further, it 
is known that an in-plane magnetic field has a strong 
effect on magnetoresistance, including hysteretic effects 
which paa be interpreted as field-induced rotation of the 
stripestH2l. Remarkably, an in-plane field also has. 



strong 

effect on the LTT-LTO transition temperaturdHi, stabi- 
lizing the LTT phase at higher T. This can be qualita- 
tively understood, if the LTT phase and charge stripes 
are intimately related: the fluctuating (LTT) stripes are 
present in the LTO phase, but the strong fluctuations 
produce an average LTO structural order; the field, by 
aligning the stripes, reduces transverse fluctuations and 
reveals the underlying LTT order. 
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FIG. 15. (a) Phase diagram of Lai.6- a; Ndo.4oSr a ;Cu04 [ 
116. Solid lines = structural transitions, with various phases 
labelled (LTO-LTT = mixed phase, LTLO = low temperature 
less orth orho mbic phase). Long dashed lines = NQR knock- 
out line [ 117 1, indicating onset of slow stripe fluctuations; cir- 
cles = charge order; dot-dashed line = long-range (incommen- 
surate) antiferromagnetic (LRAF) order; dotted line = super- 
conductivity, (b) Interpretation, in terms of dominant and 
subdominant (in parentheses) flu ctuations, (c) Spin-charge 



separation phase diagram 
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below the dashed line 



spinons are paired, below the solid line holons are bose con- 
densed. 

In summary, the data suggest the following scenario. 
The magnetic stripes are associated with a local LTO- 
type tilt distortion, the charged stripes with an LTT tilt- 
ing. At half filling, there is a simple long-range LTO or- 
der, which may be stabilized by magnetoelastic coupling 
which aligns the spins with the LTO tilt direction. The 
LTT tilting, stabilized by coupling to the VHS, is fluc- 
tuating unless pinned by ionic disorder off of the Cu02 
planes. When the LTT phase is pinned, the stripes de- 
velop long-range charge (and ultimately spin) order. If 
the stripe pinning is strong enough, the superconducting 
order can hp suppressed, predominantly due to the mag- 
netic orderlH3. Near enough to half filling, the LTO tilts 
prevail, the stripes rotate by 45°, and superconductivity 
is destroyed. 

Thus there is a three-way competition, between mag- 
netic, structural, and superconducting order. Such a 
model could explain why superconductivity is not found 
in the nickelates and manganites, where electron-phonon 



coupling is stronger. The competition between super- 
conductivity and LTT|-Qi[d^r takes place predominantly 
on the charged stripe.' 



23, and is strongly influenced 
by interlayer strains, which caut-be varied by rare earth 



substitution. Buchner, et al.crfl find a critical tilt an- 
gle $ c , such that when the average LTT tilt exceeds $ c , 
'bulk' superconductivity is absent (the Meissner fraction 
is greatly reduced), while for tilts less than <!> c supercon- 
ductivity is hardly affected by the presence of LTT order. 
Experimentally, it is often found that when LTT fluctua- 
tions are weak (e.g., in sotne compositions of LSCO and 
oxygen-doped La2Cu04llH3), magnetic order develops pt. 
actly at the superconducting transition. Julicn, et al.E^i 
have suggested that magnetic fluctuations are present at 
much higher temperatures, and that the sudden slowing 
of the fluctuations at T c is coincidental. We suggest a 
slightly different possibility: that when LTT pinning is 
particularly weak, it takes superconducting order to pro- 
vide the stripe phase with structural rigidity, and that it 
is the stiffening of the charge stripes which is reflected in 
the slowing of the magnetic fluctuations. This provides 
one more demonstration that the charged stripes can be 
either superconducting or LTT ordered. 



2. Flux Phase on Stripes 

Early studies of the flux phase generally concentrated 
on the competition between the flux and antiferromag- 
netic phases at half filling. The present analysis strongly 
suggests that the flux phase would be present instead 
in the hole doped region - that is, on the charged 
stripes. This follows both because instability is enhanced 
near a VHS, and because the correlated hopping scales 
with x; moreover, in a stripe picture, features near the 
Fermi Jewel are generally associated with the charged 
stripesEB. This coiaclusion is consistent with the experi- 
mental observationL3 that the possible flux related mag- 
netization is stronger in YBCOe+irJQr y=0.6 than for 
y=0.35 (see also Chakravarty, et aLcH3). It alsp provides 



a natural resolution of the-iee-Wen paradoxll£j discussed 



by Orenstein and MillisEHZI. Lee and Wen showed that 



they could explain the Uemura relation for underdoped 
cuprates, as long as the flux phase dispersion is inde- 
pendent of doping. The paradox is that many strong- 
coupling models expect the dispersion to renormalize to 
zero near the Mott insulator at half filling. In a stripe 
picture, this renormalization is taken as indicating that 
the fraction of material associated with charged stripes 
renormalizes to zero at half filling, whereas the. dispersion 



on a single stripe is less sensitive to dopingc 

At very low doping, flux-stripes will tend to break up 
into flux-polarons, confined to a single plaquette. These 
flux-polarons bear a close, resemblance to the skyrmions 
introduced by Gooding 



J. 2 s 



Indeed, his electronic states 
are just linear combinations of the -E„-symmetry plaque- 
tte states of Eq. ^. The difference is that his states are 
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localized around a Sr impurity, while ours are the equilib- 
rium conducting state on the hole-doped stripes. Clearly, 
at low temperatures there will be a tendency for charged 
stripes and polarons to be pinned on the Sr, greatly en- 
hancing the similarity. However, the present model more 
naturally explains stripes at higher doping levels, and the 
uniform chared phase near xq. It should be noted that 
Haskel, et al.LJ find enhanced local structural distortions 
near Sr impurities, associated with hole localizationE£j. 

If the flux phase is on the charged stripes, then hole 
doping will shift the Fermi level away from the coni- 
cal points, producing not a node but a hole pocket, as 
in Fig. [ll]). It is possible that the deviations from d- 
wave symmetry found in the underdoped regimaH3 are 
associated with this hole pocket. The apparent d-wave 
gap might then be a localization gap at the hole pocket 
Fermi surfac e, lead ing to the localization effects observed 
in resistivityi 130 lr 31 . 



1987, AndersonE^a noted that spinon and holon conden- 
sation should couple to structural distortions. Super- 
conducting order seems to be more characteristic of the 
holon condensed domain (charged stripes) than of the 
stripe phase regime. 
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3. Spinons and Holons and Stripes 

Overall, the idea that stripes involve competition 
between antiferromagnetism, d-wave superconductivity, 
and the flux phase is very attractive, and arises natu- 
rally f-ram an SO (6) symmetry (a subgroup of the SO (8) 
groupE3), in which all three states comprise a single 6- 
dimensional order parameter, Fig. [ll]. In fact, in the 
repulsive Hubbard model at half filling, when t' = 0, 
there is an extra electron-hole symmetry, leading to a 
degeneracy betyeen the flux phase and d-wave supercon- 
ductivity (dSC)E23. This is equivalent to the degeneracy 
between a CDW and s-wave superconductivity (sSC) in 
the attractive Hubbard model, since the flux phase is a 
form of d-wave CDW. In fact, the r operator of Eq. |2| 
transforms the OAF and dSC operators into the CDW 
and sSC operatorsE3. This degeneta&y is the basis for the 
SU(2) model of superconductivity^! The SO (6) model 
thus encompasses the SO(5) and SU(2) models, Fig. [n| 

There is a close similarity between the experimental 
doping dependence o f L S CQ , Fig. Jl5|a,b and the spinon- 
holon phase diagramliME 3 ^, Fig. n3c, with the spinon- 
condensed phase corresponding to the magnetic stripe 
dominated regime and the holon-condensed phase corre- 
sponding to the charged stripe dominated regime, as in- 
dicated by the experimental LTT phase, Fig. |i"5|a! Strik- 
ingly, the holon condensed regime is where the holes have 



a flux-phase Fermi surfaceEM Indeed, it would seem 



that spin-charge separation can be interpreted as a the- 
ory of stripe phase formation. Thus, the spinon pair- 
ing involves the exclusion of spin deviations from an an- 
tiferromagnetic domain onto an antiphase domain wall, 
while holon condensation involves the ordering which sta- 
bilizes the charged domains (the flux phase). The spin 
gap would then correspond to a quantum confinement 
energyfcia. When both spinons and holons are condensed 
together, well-defined stripe order results. As early as 



FIG. 16. Six dimensional order parameter in SO(6) theory 
with three components of magnetic order (AFM), two of 
d-wave superconductivity (dSC), and one of flux phase (flux). 
In a stripe picture, the shaded plane would correspond to the 
charged stripes. 
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